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ABSTRACT 

We describe the data processing pipeline employed by the Low Frequency Instrument (LFI) data processing centre (DPC) to create and characterize 
the full-sky maps based on the first 15.5 months of operations at 30, 44 and 70 GHz. In particular, we discuss the various steps involved in reducing 
the data, starting from telemetry (TM) packets through to the production of cleaned calibrated timelines and calibrated frequency maps. Data 
are continuously calibrated using the modulation induced on the mean temperature of the cosmic microwave background radiation by the proper 
motion of the spacecraft. Sky signals other than the dipole are removed by an iterative procedure based on simultaneous fitting of calibration 
parameters and sky maps. Noise properties are estimated from time-ordered data after the sky signal has been removed, using a generalized least 
square map-making algorithm. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, 
minimizing the variance of correlated noise. Noise covariance matrices, required to compute statistical uncertainties on LFI and Planck products, 
are also produced. Main beams are estimated down to the a -20 dB level using Jupiter transits, which are also used for the geometrical calibration 
of the focal plane. 

Key words. methods:data analysis; cosmology: cosmic microwave background 



1. Introduction 

This paper, one of a set associated with the 2013 release of 
data from the Planck 1 mission (Planck Collaboration I 2013), 



* Corresponding author: A. Zacchei, zacchei@oats.inaf.it 
1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 



describes the LFI (low frequency instrument) data processing 
which supports the first Planck cosmological release based on 
the nominal Planck survey (15.5 months of observation). This 
paper represents an updated version of the LFI data process- 
ing description (Zacchei et al. 2011) which was part of the 
first wave of early astrophysical results published in early 201 1 
(Planck Collaboration, 2011h-z). This work describes the over- 



telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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all data flow of the pipeline implemented at the LFI DPC, from 
instrument scientific telemetry and housekeeping data to fre- 
quency maps, as well as the test plan applied to validate the 
data products. Due to the complexity of the analysis process, 
three companion papers and one explanatory document will pro- 
vide further details on specific critical aspects of the data anal- 
ysis and products delivered, and they will be referenced here 
as appropriate. In Planck Collaboration III (2013) we discuss 
systematic effects in detail and present an overall error budget. 
In Planck Collaboration V (2013) we describe photometric cal- 
ibration, and provide details of the methods implemented and 
the related uncertainties. Next the determination of the detector 
main beams and their uncertainties, based on in-flight planet- 
crossing measurements, is presented in Planck Collaboration IV 
(2013). Finally in Planck Collaboration ES (2013) description 
of the products delivered is detailed. When needed, the justifi- 
cation for the choice of the algorithm used in the pipeline will 
be discussed in the relevant companion papers. All the main in- 
formation and reference tables discussed in detail in companion 
papers (e.g., noise characterization, beam parameters, etc.) will, 
however, be reported and summarized in this work. 

2. LFI in-flight behaviour and operations 

The Planck LFI instrument is described in detail in Bersanelli 
et al. (2010) and Mennella et al. (2010). It comprises eleven 
Radiometer Chain Assemblies (RCA's), two at 30 GHz, three 
at 44 and six at 70 GHz, each composed of two independent 
pseudo-correlation radiometers sensitive to orthogonal linear po- 
larizations. Each radiometer has two independent square-law 
diodes for detection, integration and conversion from RF signal 
into DC voltages. The LFI focal plane is cryogenically cooled to 
20 K, while the pseudo-correlation design uses internal, black- 
body, reference loads cooled to ~ 4.5 K (4KRL). The radiome- 
ter timelines are produced by taking differences between the 
signals from the sky, V s k y , and from the reference loads, V re f. 
Radiometer balance is optimized by introducing a gain modu- 
lation factor, typically stable within 0.04% throughout the mis- 
sion, which greatly reduces 1 // noise and improves immunity 
to a wide class of systematic effects (see Mennella et al. 2011, 
for further details). During the entire nominal survey, the behav- 
ior of all 22 LFI radiometers was highly stable, with 1 // knee 
frequencies of order 50 mHz and white noise levels unchanging 
within 0.5%. Because of the stability of the instrument, only a 
very limited number of corrections were required in the LFI data 
processing, either at the TOI (time ordered information) and map 
level. 

2.1. Operations 

During the period of observation the entire satellite was very sta- 
ble (Planck Collaboration 1 2013), with a single exception. Three 
months before the end of the nominal mission it was necessary, 
as scheduled, to switch from the nominal to the redundant sorp- 
tion cooler. This operation, briefly described below, was visible 
in the LFI scientific data, but the effect on the temperature power 
spectrum was estimated to be negligible (see Sect. 2.2) 

2.2. Switch-over from nominal to redundant Sorption Cooler 

The 20 K cooling on Planck is provided by the sorption cooler 
system (SCS). This cooler uses six metal hydride compressor 
elements to produce high-pressure hydrogen that is circulated 



through a Joule-Thomson valve to provide 1 W of cooling at 
20 K. Gas compression is achieved by heating a single compres- 
sor element to 440 K and a pressure of 30 Bar. After expansion 
through the Joule-Thomson valve, the gas is recovered by three 
compressor elements at 270 K and 0.3 Bar. To reduce power con- 
sumption, gas-gap heat switches are used to isolate the com- 
pressor elements from the radiator during gas compression. Two 
sorption coolers were flown on the Planck spacecraft to meet 
mission lifetime requirements. A switchover procedure was de- 
veloped to change between the operating cooler and the redun- 
dant cooler. In early August of 2010, one of the gas-gap heat 
switches for a compressor element failed on the active cooler. 
Although the SCS can operate with as few as four compressor 
elements, it was decided to implement the switchover proce- 
dure and activate the redundant cooler. On 11 August at 17:30 
GMT, the working cooler was commanded off, and the redun- 
dant cooler was switched on. Adequate cryogenic cooling was 
restored in about 1 hour. Return to thermal stability, however, 
took 48 hours. After thermal stability of the SCS was restored, 
anomalous temperature fluctuations were observed on the LFI 
focal plane. These excess fluctuations are thought to be due to 
sloshing of liquid hydrogen remaining in the cooler that had 
been switched off. This was due to the loss of storage capac- 
ity of the metal hydride of the inactive cooler, which prevented a 
reduction of the pressure to a level that prevented liquid forma- 
tion. While these fluctuations produced a measurable effect in 
the LFI data, their propagation to the temperature power spec- 
trum is more than two orders of magnitude below the CMB sig- 
nal (Planck Collaboration III 2013). Furthermore, after the end 
of the nominal mission in February 201 1, these fluctuations were 
reduced to a much lower level. More details of these issues will 
be discussed in a future paper. 

2.3. Instrument performance update 

Table 1 gives a top-level summary of instrument performance 
parameters measured in flight during the nominal data period. 
Optical properties have been successfully reconstructed using 
Jupiter transits (Planck Collaboration IV 2013) and the main pa- 
rameters are in agreement with pre-launch and early estimates 
(Mennella et al. 201 1). The white noise sensitivity and parame- 
ters describing the 1 // noise component are in line with ground 
measurements (Mennella et al. 2010) and agree with the values 
published in the early paper (Mennella et al. 201 1) Photometric 
calibration based on the CMB dipole, yields an overall statistical 
uncertainty of 0.25% (Planck Collaboration V 2013). Variations 
due to slow instrumental variations are traced by the calibra- 
tion pipeline, yielding an overall uncertainty ranging between 
0.1% and 0.2%. The residual systematic uncertainty is linked to 
the frequency and varies between 21 muK and 6 /jKqmb (Planck 
Collaboration III 2013). 

3. LFI Data processing overview 

The processing of LFI data can be logically divided into a se- 
ries of levels shown schematically in Fig. 1. The process starts 
with the Level 1 (LI); this has the responsibility to retrieve all 
necessary information from packets and auxiliary data received 
each day from the Mission Operation Center and to transform 
the scientific packets and housekeeping data into a form that is 
manageable by the Level 2 (L2) scientific pipeline. The Level 2 
pipeline queries the database where the LI data has been stored 
and, using scientific and housekeeping information, does the fol- 
lowing 
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Table 1. Summary of main LFI performance parameters. 



Parameter 


j\J OrlZ 


A A ("'11/ 


/U OrlZ 


Centre frequency [GHz] 


28.4 


44.1 


70.4 


Scanning Beam " FWHM ['] 


33.16 


28.09 


13.08 


Scanning Beam " Ellipticity 


1.37 


1.25 


1.27 


Effective Beam * FWHM ['] 


32.34 


27.12 


13.31 


White noise sensitivity (map level) c [/jKcmb] 


9.2 


12.5 


23.2 


White noise sensitivity (from timelines) [yiiKcMB s ] 


148.5 


173.2 


151.9 


/knee [mHz] 


114.5 


45.7 


20.2 


1// slope 


-0.92 


-0.90 


-1.13 


Overall calibration uncertainty [%] 


0.82 


0.55 


0.62 


Systematic effects uncertainty [/jK C mb] 


21.02 


5.61 


7.87 



11 FWHM and ellipticity from scanning beam evaluated fitting Jupiter directly in the timelines. 

b FWHM from effective beam estimated from the main beam solid angle of the effective beam, see Sec. 6.2. Those are the values used in the 

source extraction pipeline (Planck Collaboration XXVIII 2013). 
c White noise per pixel computed from half-ring difference maps. Those values are within the 1% with the white noise sensitivity computed 

directly on the timelines, taking in account the actual mission time minus the manoeuvres and bad-science flagged data. 



- build the LFI RIMO (reduced instrument model) that con- 
tains all the main characteristics of the instrument; 

- remove AD/C non linearity and 1 Hz spikes at diode level 
(see Sects. 4.2 and 4.3); 

- compute and apply the gain modulation factor to minimize 
the 1// noise (see Sect. 4.4); 

- combine the signals from the diodes (see Sect. 4.5); 

- compute the corresponding detector pointing for each sam- 
ple, based on auxiliary data and beam information (see 
Sect. 5); 

- calibrate the scientific timelines to physical units (Kcmb) fit- 
ting the dipole convolved with the 47r beam representation 
(see Sect. 7); 

- remove from the scientific calibrated timeline the dipole con- 
volved with the 47r beam representation; 

- combine the calibrated TOIs into aggregate products such as 
maps at each frequency (see Sect. 9). 

The Level 3 collects instrument-specific (from both HFI 
and LFI) L2 outputs and derives various products. For in- 
stance, component-separation algorithms transform frequency 
maps into maps of separated astrophysical components; source 
detection algorithms create catalogues of different classes of ob- 
ject and sources; and finally a likelihood code assesses the match 
between a cosmological and astrophysical model, on the one 
hand, and the frequency maps on the other. 

4. Time ordered information (TOI) processing 

In the LI pipeline the telemetry data are received at the DPC 
(data processing center) as a stream of packets which are handled 
automatically by different steps in the Level 1 : 

- uncompress the retrieved packets; 

- de-quantize and de-mix the uncompressed packets to retrieve 
the original signal in ADU (analog-to-digital units); 

- using a conversion factor stored in the packet header to trans- 
form ADU data in to volts; 

- cross correlate time information to univocally time stamp 
each sample; 

- store the resulting timelines into a database interface to the 
Level 2 pipeline. 



We made no change in Level 1 software during the mission. 
Detailed information on how each of the steps listed above were 
applied is provided in the early paper by Zacchei et al. (201 1). In 
order to avoid strong gradients in the signal and signals that do 
not project correctly in the maps, we establish a procedure able 
to flag a single scientific sample; this is described in Sect. 4. 1 . 

4. 1 . Input flags 

For each sample we define a 32-bit flag mask to identify potential 
inconsistencies in the data and to enable the pipeline to skip that 
sample or handle it differently. The TOI from all LFI detectors is 
archived in the Level 1 database, and regularly checked to iden- 
tify and flag events that can affect the scientific analysis. These 
events include missing or anomalous data, and data acquired dur- 
ing the manoeuvres regularly performed to repoint the telescope 
according to the Planck scanning strategy. Table 2 summarizes, 
for the period of the nominal mission, the percentage of time 
associated with these events in the LFI TOIs. The table also re- 
ports the total percentage of Level 1 TOIs usable in the scientific 
analysis. Most of the missing data results from telemetry packets 
where the arithmetic compression performed by the SPU is in- 
correct, causing a decompression error. They are very rare, and 
hence they cause a completely negligible impact on the scientific 
analysis. For instance, for the entire 70 GHz frequency channel, 
the total amount of missing data corresponds to 130 lost seconds 
in 15 months. The instrument team performs a daily check of 
the LFI data retrieved during the daily tele-communication pe- 
riod with the satellite; the data covers an entire operational day. 
Part of this analysis consists in identifying, for each LFI detector, 
time windows where either the total power signal or the differen- 
tiated signal shows anomalous fluctuations or jumps. Depending 
on the characteristics of the anomaly identified, a time window 
can be flagged as unusable for science. Currently, the criteria 
defined to flag time windows as unusable include: 

- time windows where a gain change in the data acquisition 
electronics caused a saturation of the sky or reference load 
signals; 

- abrupt changes in voltage output with slow recovery (>1 m), 
caused by gain fluctuations in the back-end module ampli- 
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Raw data from LI 









AD/C& 




LI raw TOI 


1Hz spikes 




diode 



AD/C Non linearity correction 
& 

1 Hz Spikes Removal 



Check, Fill Gaps & flag 



Subsampled 
data 



Sky/load mean & GMF evaluation 



/ diode 7 
/ weights / 




LEVEL 2 pipeline 

Fig. 1. Schematic representation of Level 1, Level 2 and pointing pipeline of LFI DPC 



fier, induced by electrical or thermal variations, which also 
generate discontinuities in the differentiated signal; 

- short, abrupt changes in voltage output caused by fluctua- 
tions in the low noise amplifiers in the front end module, 
which produce asymmetries between the sky and reference 
load signals and possibly first order effects in the differenti- 
ated signal; 

- permanent changes in the output voltage caused by a perma- 
nent change at the front end module (amplifier bias or focal 
plane unit temperature) or back end module (temperature or 
HEMT gain variations); in such cases, only a small time win- 
dow around the discontinuity is flagged as unusable; 

- short time windows (< 1 m) where the total power signal 
shows pop-corn noise on one or both detectors, due to vari- 
ations in the bed end diode or in the front end low noise 
amplifiers 

In Table 2, the row labelled as "Anomalies" reports the per- 
centage of observation time flagged as unusable for these rea- 
sons in the scientific analysis. The almost 1.0% shown for the 
44 GHz channel corresponds to a total time of 1 1 3 hours. Finally, 
the times of manoeuvres and stable pointing periods are recov- 
ered from the attitude history files provided by the Planck flight 
dynamics team. Detector samples corresponding to manoeuvres 
are flagged out, in order to skip them during the present analysis. 



Table 2. Summary of the percentage of the LFI observation time 
spent in manoeuvres, lost due to missing or unusable data and ef- 
fectively used in the scientific analysis, calculated by frequency 
channel. 





30 GHz 


44 GHz 


70 GHz 


Missing [%] 


0.00014 


0.00023 


0.00032 


Anomalies [%] . . . . 


0.82220 


0.99763 


0.82925 


Maneuvers [%] . . . . 


8.07022 


8.07022 


8.07022 


Usable [%] 


91.10744 


90.93191 


91.10021 



Tasks within the Level 2 pipelines both fill gaps in the data 
with artificial noise and flag them properly. Other tasks locate 
planet transits and moving objects within the Solar system, again 
flagging samples affected by such observations. 

4.2. ADC linearity correction 

The analogue to digital converters are the part of the data acqui- 
sition electronics responsible for producing the digital presenta- 
tions of the analogue detector voltages which are then processed 
on-board by the radiometer electronics box assembly. Since they 
are directly involved with the signal power, their linearity is as 
important as that of the receivers and detectors, with any de- 
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parture appearing as a distortion in the system power response 
curve. In differential measurements such as those carried out by 
the Planck LFI instrument, very small localized glitches in this 
curve can have a large impact, since the calibration factor de- 
pends on the gradient of the response curve at the point at which 
the differential measurements are made. This effect is described 
in detail in the systematics paper Planck Collaboration III (2013) 
and the impact on calibration in Planck Collaboration V (2013). 

The effect is observed in some LFI radiometer data, appear- 
ing as gain variations seen at particular detector voltages. This is 
shown for the worst affected channel RCA2501, in Fig. 2, where 
the upper plot shows the measured voltages of the sky and ref- 
erence loads (in red and blue respectively) and the lower plot 
the percentage variations of gain (green) and noise on sky and 
reference voltages (red and blue). The range of the upper plot is 
matched to that of the lower plot and so for normal gain vari- 
ations the same pattern should be seen for both. That is clearly 
not the case. When the sky signal falls close to 0. 186 V, as shown 
by the horizontal dashed lines, both the inverse gain and the sky 
"white noise" estimates show anomalies in the same time in- 
terval (the time interval affected is indicated by vertical dashed 
lines in the Fig. 2). The same anomalous behavior of the ref- 
erence "white noise" signal and inverse gain is seen in two in- 
tervals when the reference signal falls close to either 0.197 V 
or 0.202 V. Outside of these limited ranges, the variations in all 
plotted signals track one another, such as the feature at 192 days 
after launch in the sky voltage or the drop at day 257 when the 
transponder was turned permanently on. 

The response curves can be reconstructed by tracking how 
the noise amplitude varies with the apparent detector voltage in 
the TOI. The radiometers are assumed to be stable and the in- 
trinsic thermal noise can be taken to be constant in terms of tem- 
perature, so any voltage variations are then assumed to be due to 
both gain drift and ADC effects. The method for this correction 
is set out in appendix A of Planck Collaboration III (2013). 




100 150 200 250 

Days after launch 



300 



Fig. 2. Effect of ADC non-linearities on time ordered data of 
one of the 44 GHz diode. The upper plot shows the recorded 
detector voltage (sky in red and reference in blue). The hori- 
zontal shaded bars show the voltage range affected by the ADC 
non-linearities and the vertical bars the time range affected. The 
lower plot shows the percentage variation of the inverse of the 
gain factor from the dipole gain (green) and the "white noise" 
estimates on the sky and reference voltages (again sky red and 
reference blue). The gain estimates have been smoothed by a 
three day moving mean and the noise by one day. 



ing pipeline, is to take the ratio of DC levels from sky and 
load outputs obtained by averaging the two time streams, i.e., 
R = <Vsky)/<Vioad). Then 



4.3. Corrections for electronic spikes 

This effect is caused by interaction between the housekeeping 
electronics clock and the scientific data line in the on-board data 
acquisition system. The spikes are synchronous with the on- 
board time with no changes in phase over the entire acquisition 
period, allowing the construction of dedicated templates that are 
then removed from the timelines. Frequency spikes are present 
in the output at all frequency but affect mostly the 44 GHz due to 
the high gain in these detectors: consequently, electronic spikes 
are removed only in this channel. This process and the evaluation 
of the effect at map level is described in Planck Collaboration III 
(2013). 

4.4. Demodulation: Gain modulation factor estimation and 
application 

Each diode switches at 4096 Hz (Mennella et al. 2010) between 
the sky, V s k y , and the 4K reference load, Vi oa d- Vsky and Vi oa d 
are dominated by 1 // noise, with knee frequencies of tens of 
hertz. This noise is highly correlated between the two streams, a 
result of the pseudo-correlation design (Bersanelli et al. 2010), 
and differencing the streams results in a dramatic reduction of 
the 1/f noise. To force the mean of the difference to zero, the 
load signal is multiplied by the GMF (gain modulation factor), R, 
which can be computed in several ways (Mennella et al. 2003). 
The simplest method, and the one implemented in the process- 



AV(f) = V sky (0 - — ^" Vload(0 . (i) 
{ Vioad) 

The R value was computed from unflagged data for each point- 
ing period and then applied to create the differenced timelines. 
The R factor has been stable over the mission so far, with overall 
variations of 0.03-0.04%. Full discussion regarding the theory 
of this value is reported in Mennella et al. (201 1) 

4.5. Combining diodes 

The receiver architecture is symmetric, with two complemen- 
tary detector diodes providing output for each receiver chan- 
nel. As described in Seiffert et al. (2002) and Mennella et al. 
(2010), imperfect matching of components limits the isolation 
between the complementary diodes of the receivers to between 
-10 and -15 dB. This imperfect isolation leads to a small anti- 
correlated component in the white noise. We perform a weighted 
average of the time ordered data from the two diodes of each 
receiver just before the differentiation. This avoids the compli- 
cation of tracking the anticorrelated white noise throughout the 
subsequent analysis. We treat the combined diode data as the raw 
data, and calibration, noise estimation, mapmaking etc. are per- 
formed on these combined data. We use inverse noise weights 
determined from an initial estimate of the calibrated noise for 
each detector. The weights, reported in Table 3, are kept fixed 
for the entire mission. 
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Table 3. Diode weights applied (first three decimals), perfect 
instrument would have 0.500 weights for both diodes. 



M-00 M-01 S-10 S-ll 



LFI 18 0.567 0.433 0.387 0.613 

LFI 19 0.502 0.498 0.551 0.449 

LFI20 0.523 0.477 0.477 0.523 

LFI21 0.500 0.500 0.564 0.436 

LFI22 0.536 0.464 0.554 0.446 

LFI23 0.508 0.492 0.362 0.638 

LFI24 0.602 0.398 0.456 0.544 

LFI25 0.482 0.518 0.370 0.630 

LFI26 0.593 0.407 0.424 0.576 

LFI27 0.520 0.480 0.485 0.515 

LFI28 0.553 0.447 0.468 0.532 



5. Pointing 

Proper pointing reconstruction is critical and has a direct impact 
in the determination of an accurate photometric calibration. The 
pointing for each radiometer P a d(i) at time t is given by 

PyadiO = ^?Ecl,Body(0^?Body,radCz- (2) 

The ^Body.rad matrix encodes the orientation of the beam pat- 
tern with respect to the body reference frame, i.e., the reference 
frame defined by the spacecraft structure. We adopt the conven- 
tion that in the reference frame of the beam the optical axis is 
aligned with e z . The !??Body,rad is parameterized by a set of ro- 
tation angles in the RIMO derived from flight data and ground 
based measurements. The 9?Eci,Body(0 is derived by time inter- 
polation of quaternions distributed in the attitude history files. 
The spacecraft attitude is determined by the Planck Star Tracker 
and during periods of stability between maneuvers is sampled at 
8 Hz, much lower than the LFI sampling frequency. Despite its 
simplicity, this formula synthesizes a large amount of informa- 
tion on the satellite and a long chain of transformations between 
reference frames, each one being a possible source of systematic 
errors. Indeed, even a small aberration compared to the beam 
size can introduce significant photometric effects if the gradient 
of the temperature field is large enough. The two most impor- 
tant sources of aberration identified and corrected in PlanckfLFl 
are stellar aberration and the apparent change in wobble angles 
likely produced by thermal deformations of the star tracker sup- 
port. 

5. 1 . Stellar Aberration 

The star tracker system provides reconstructed astrometric atti- 
tude for the Planck spacecraft in the Solar System baricentric 
reference frame (SSB). However the effective pointing direction 
is aberrated by the well known effect of orbital stellar aberration 
due to the combination of Planck orbital motion and finite speed 
of light. In the non-relativistic case stellar aberration follows the 
well known equation 

P' = {P + Vpi anck /c)/ \P + Vpi anck /c| , (3) 

where P' is the aberrated pointing direction, Vpi anc k is the orbital 
velocity of Planck in the SSB reference frame, and c the speed of 
light. From this formula the deflection angle dp = arccos P ■ P 
can be derived. Since Planck moves at about 30Kms~' in the 
Solar System, its motion is nearly coplanar to the ecliptic and 



scan circles are nearly normal to it: 5p < 20.6", the worst de- 
flection occurs near the ecliptic poles. If left uncorrected, this 
aberration would distort the maps and it will be more similar 
to a seasonal shift near the equator and to a "blurring" near 
the ecliptic poles. Indeed a set of accurate simulations shows 
that the distortion radius is maximal at the ecliptic poles, i.e., 
at mid Galactic latitudes (Z G ai,&Gai) = (96.384°, 29.81 1°) and 
(faai.&Gai) = (276.384°, -29.8 11°), and that it decreases toward 
the ecliptic down to a minimum of about 0. 1 ° on the ecliptic. The 
boundary of the worst affected region, i.e., the region in which 
the distortion radius is at least half the polar value, is roughly 
a ring centred on the poles with radius about 60°. There are 
some variations in the radius and in the longitudinal shape of 
the boundary, both smaller than a few degrees, due to the scan- 
ning strategy and also to the different angular distances from the 
spin axis of the various feed-horns. 

5.2. Wobble Angles 

Wobble angles describe the unavoidable misalignment of the 
body reference frame with respect to the reference frame de- 
fined by the satellite principal inertial axis. The nominal spin 
axis for the satellite is nearly a half degree away from the princi- 
pal moment of inertia and the effective scan circles are about half 
a degree smaller than the nominal ones (Planck Collaboration I 
2011). Wobble angles and their variations in time, either real 
or apparent, are measured by careful modelling of the observed 
Planck attitude dynamics included in the attitude history files. 
Planck Collaboration I (201 1) reported an apparent variation of 
the wobble angles likely produced by thermoelastic deforma- 
tions which change the relative orientation of the star tracker 
with respect to the body reference frame. The change was de- 
tected in scans of Jupiter. Since this variation is rigidly trans- 
ported by the rotations of spacecraft body, its effect will be 
largely averaged out near the poles and will be maximal near 
the ecliptic, the opposite of the stellar aberration effect. Of the 
three angles that describe the wobble, 4>\ has largely negligible 
effects and ^3 is badly determined, so the LFI pipeline applies 
a correction for only the 1^2 variation. Its effect is like an appar- 
ent change of the angular distance between the telescope and the 
spin axis so that a spurious change of this angle will be equiva- 
lent to an apparent change of scan circles radii of ±0.1', giving 
an equivalent displacement in pointing between consecutive sur- 
veys of 0.2'. 

6. Main Beams and the geometrical calibration of 
the Focal Plane 

Through measurements of Jupiter crossing we are able to mea- 
sure the beam profiles (Planck Collaboration IV 2013) and di- 
rections. We use this information to determine the focal plane 
geometry, which is the pointing direction and its related polar- 
ization for each beam centre in the field of view. The origin of the 
focal plane is taken to coincide with the optical line of sight de- 
fined in Tauber et al. (2010). In the LFI RIMO the centre of beam 
is expressed by four numbers, namely 6 UV ; <p uv ; i/r ul ,; i// po i; see 
Planck Collaboration ES (2013) for the definition of these an- 
gles. Only 8 av and <p uv , which are the beam pointing in spherical 
coordinates referred to the line of sight, can be determined with 
Jupiter observations; the polarization orientation of the beams, 
defined by if/ m , + if> po i, is not estimated from flight data but is 
derived from main beam simulations based on ground measure- 
ments. 
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Table 4. Approximate dates of the Jupiter observations. 



Jupiter crossing Date OD 



Scan Jl 21/10/09 - 05/11/09 161 - 176 

Scan J2 27/06/10 - 12/07/10 410 - 425 

Scan J3 03/12/10 - 18/12/10 569 - 584 

Scan J4 30/07/ 1 1 - 08/08/ 11 808-817 



The calibration of the beams and their locations was accom- 
plished using the four Jupiter crossings which occurred in the pe- 
riods reported in Table 4, following the same procedure already 
described in Zacchei et al. (201 1) and Mennella et al. (201 1). 

For each beam, the pointing is determined by the location of 
the maximum of the elliptical Gaussian fit to that beams. This 
was done for each beam in each single scan whose results are 
reported, with errors, in Planck Collaboration IV (2013). 

In addition, the beams were stacked in pairs (J1J2 and J3J4) 
and all together (J1J2J3J4) in order to improve the signal to noise 
ratio of the measurements. Before the stacking each beam was 
artificially repointed along the direction given by the arithmetic 
average of the centre of each single beam to be stacked. Then a fit 
was performed again on the stacked beams and the resulting pa- 
rameters recorded. For single scans it has been found that there is 
an agreement within 2"in the pointing direction between Jl and 
J2. The same agreement occurs between J3 and J4. In contrast, 
a ~ 15 "systematic deviation of the beam centre was detected 
when comparing Jl J2 to J3J4. Figure 3 reports the reconstructed 
beam positions and errors in the line of sight frame both mag- 
nified by a factor of 100. The shift is evident for the 70 GHz 
beams as well as in all the J1J2 and J3J4 stacked beam centres. 
It became manifest when we stacked the beams to improve the 
determination of their centres. The change in the location has 
been found mainly in the scan direction (i.e., v-coordinate). To 
account for this pointing shift we apply two pointing solutions 
for LFI. The first focal plane calibration is valid from OD91 to 
OD540 and it is based on the Jl J2 beam pointing determination. 
The second calibration is valid from OD541 to OD563 and it is 
based on the J3J4 beam pointing calibration. The reconstructed 
angles are reported in Table 5 

6.1. Scanning beams 

Scanning beams are defined as the measured beams in flight. 
They include the smearing effect due to the satellite spinning, 
and the effects of the radiometers, such as the band shape. With 
four Jupiter transit measurements we were able to reconstruct the 
beam shape down to -20 dB from peak at 30 and 44 GHz, and 
down to -25 dB at 70 GHz. From the beam shape we estimated 
the main beam parameters using a bivariate Gaussian fit on the 
four stacked beams (J1J2J3J4). The fitting procedure, described 
in Zacchei et al. (201 1), was slightly modified to correct for off- 
set in the data and to avoid noise contamination. We refer to the 
companion paper on LFI beams (Planck Collaboration IV 2013) 
for details on procedures and results. Here we report in Table 6 
the average values of the FWHM and ellipticity with errors. 

6.2. Effective Beams 

The effective beam is the average of all scanning beams pointing 
at a certain direction within a given pixel of the sky map with 
the Planck scanning strategy applied. We compute the effective 



beam at each LFI frequency scanning beam and scan history in 
real space using the FEBeCoP (Mitra et al. 201 1) method. Details 
of the application of FEBeCoP to Planck data will be discussed in 
a future paper. Effective beams were used to calculate the effec- 
tive beam window function as reported in Planck Collaboration 
IV (2013) and in the source detection pipeline necessary to gen- 
erates the PCCS catalogue (Planck Collaboration XXVIII 2013) 
We list in Table 7 the main parameters computed with FEBeCoP, 
Note that the FWHM and ellipticity in Table 7 differs slightly 
from the values reported in Table 6; this results from the different 
way in which the Gaussian fit was applied. The scanning beam 
fit was determined by fitting the profile of Jupiter on timelines 
and limiting the fit to the data with a signal above the 3cr level 
from the noise, while the fit of the effective beam was computed 
on GRASP maps projected in several position of the sky (Planck 
Collaboration IV 2013); the latter are less affected by the noise. 

7. Photometric Calibration 

In this section we discuss the procedures used to perform a pho- 
tometric calibration of the PlanckfLFl radiometric data, that is, 
the conversion of the time-ordered streams of voltages into time- 
ordered streams of thermodynamic temperatures. The conver- 
sion is modelled by the following equation: 

V = G x (r sky + 

* noise ), (4) 

which relates the voltage V measured by the ADC with the 
temperature r s k y , obtained by convolving the sky temperature 
with the beam response of the instrument at a given time, and 
the noise temperature r no j se of the radiometer. In general we 
are interested in K = G , as the purpose of the calibration is 
to convert V back into a temperature. As described in Planck 

LFI M-Beams boresight direction in Planck FOV 

error boxes and positions magnified by a factor 100 
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Fig. 3. Main beam pointing directions measured with the first 
four Jupiter crossings: single scans (yellow, light red, green, 
light blue); first and second stacked scans (red); third and fourth 
stacked scans (blue); and four stacked scans (grey) are reported. 
The colored boxes refer to the measured uncertainties magnified 
by a factor of 100. The differences in pointing were normalized 
to the Jl measurements, and were magnified by the same factor 
of 100. 
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Table 5. Focal plane geometry 
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11 Beam pointing reconstructed using the first two Jupiter transits (Jl and J2). 
b Beam pointing reconstructed using the last two Jupiter transits (J3 and J4). 
c Polarization orientation of the beam measured during ground test. 



Table 6. LFI beam FWHM and ellipticity measured in flight 
from four Jupiter passes. 
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Uncertainty 13 
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Ellipticity 


Uncertainty 13 


70 GHz mean 
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1.27 




LFI18 . 
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0.03 
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0.01 
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11 The square root of the product of the major axis and minor axis 
FWHMs of the individual horn beams, averaged between M and S 
radiometers. 

b The standard deviation of the mean of the lcr statistical uncertainties 
of the fit. Although a small difference between the M and S beams 
caused by optics and receiver non-idealities can be expected. 

c Ratio of the major and minor axes of the fitted elliptical Gaussian. 

Collaboration V (2013), two different algorithms were used for 
calibrating the LFI radiometers in this data release: 

1. For the 44 and 70 GHz radiometers, we have used a tech- 
nique called OSG (optimal search of gain), which is similar 
to the one used by WMAP (Hinshaw et al. 2009). It is based 



on the fitting of the radiometric signal against the expected 
dipolar anisotropy induced by the motion of the spacecraft 
with respect to the CMB rest frame. 
2. For the 30 GHz radiometers, we have used a technique that 
combines the knowledge of the dipolar anisotropy (as above) 
then additionally takes into account the observed fluctuations 
in the measurement of the signal of the 4 K reference loads 
onboard Planck. 

The overall accuracy in the calibration is reported in Table 1 . 
The reasons why we used two different algorithms are discussed 
in the Planck Collaboration V (2013). In the following sections 
we only outline the algorithms. 

7.1. Iterative calibration 

The iterative calibration algorithm used in the Planck LFI can be 
summarized by the following points: 

1 . We combine the speed of the spacecraft with respect to the 
Sun, Vpi anc k, and the speed of the Sun with respect to the 
CMB, vsun ■ To take in account the beam shape we rotate 
the resulting vector in the beam reference frame using the 
detector pointings and estimate the direction of the dipole as 
the dot product between the rotated velocity and previously 
stored beam parameters. The dipole is then evaluated con- 
sidering the relativistic correction: 



AT = 7cmb 



1 



r(i -/?cos0) 



- 1 



(5) 



where 7cmb = 2.7255 K is the temperature of the CMB 
monopole and cos 9 is the direction computed above. We 
produce discrete time ordered data (TOD) of the expected 
overall dipole signal for each sample in a pointing period. 
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Table 7. Mean and standard deviation of FWHM, ellipticity, orientation, and solid angle of the FEBeCop effective beams computed 
with the GRASP beam fitted scanning beams. FWHM e ff is the effective FWHM estimated from the main beam solid angle of the 
effective beam, Q e ff = mean(Q). 



Frequency mean(FWHM)' cr(FWHM)' mean(e) cr(e) mean(i^)° o"(i/0° mean(fl) [arcmin 2 ] <x(Q) [arcmin 2 ] FWHM c() 



70 13.252 0.033 1.223 0.026 0.587 55.066 200.742 1.027 13.31 

44 27.005 0.552 1.034 0.033 0.059 53.767 832.946 31.774 27.12 

30 32.239 0.013 1.320 0.031 -0.304 55.349 1189.513 0.842 32.24 



2. Using pointing information, we project both V, and AT"; on a 
HEALPix map (N s id e = 256). Multiple hits on the same pixels 
are averaged in both cases. The result is a pair of maps, V™ ap 
and A7 , ™ ap , with k being the pixel index 2 . 

3. We use weighted least squares to estimate K in eq. 4 from 
the correlation between the signal in volts, Vf ky , with A7? ky : 



(6) 



where K and e are the parameters used in the fit. Each sample 
k is weighted according to the number of hits per pixel. In 
computing the fit, we use a frequency dependent mask to 
avoid those pixels where a strong non Gaussian signal other 
than the dipole is expected, i.e., point sources and the Galaxy. 

4. The main source of uncertainties in the fit using the dipole is 
given by the cosmological (CMB) signal. To improve the re- 
sult we calibrate the data using Kj and e,-, remove the dipole 
convolved with the beam and build a map, which represents 
an estimation of the cosmological signal. To reduce the effect 
of the noise we combine data streams from both radiometers 
of the same horn. Then we remove the estimated signal from 
the data and build a map, apply a simplified destriping algo- 
rithm and use the results to correct both Kj and e,-. We iterate 
the procedure until convergence. The result of this process is 
a set of gains, Kf n , and offsets, e" el . 

5. An adaptive low-pass filter based on wavelets is applied to 
the vectors Kf a and e" er in order to reduce high-frequency 
noise, particularly near the regions where the spacecraft is 
unfavorably aligned with the dipole. 

7.2. Calibration using 4K reference load signal 

To calibrate the PlanckfLFl 30 GHz radiometers, we used a dif- 
ferent calibration scheme based on the signal measuring the tem- 
perature of the 4K reference loads. This calibration has the ad- 
vantage of being less dependent on optical systematics such as 
far sidelobes (Planck Collaboration III 2013), at the expense of 
being more sensitive to systematics in the radiometers, like ADC 
non-linearities (Planck Collaboration V 2013). The algorithm 
can be summarized in the following points: 

1 . For each pointing period i, a set of gains Kf a is estimated 
using the iterative procedure described in section 7.1. 

2. The values of Kf er are used to estimate the value of the con- 
stant Kq in the equation 



Kf ei = K x 



2- 



ef \ 



V ref ,0j 



(7) 



2 Most of the pixels in the maps are not set, as during one pointing 
period the beam paints a thin circle in the sky. We assume hereafter that 
the index k runs only through the pixels which have been hit at least 



where V[ ef is the average value of the 4 K reference load sig- 
nal (in volts) over the r'-th pointing period and V** , is a 
voltage representative of the value of V™ f over the whole 

mission. (We use VJ ef = / VT ef y.) The constant Kq is esti- 
mated using a weighted, one-parameter linear least squares 
fit, where the weights are chosen to be proportional to the ex- 
pected amplitude of the dipole-like signal in the sky, A!T dlp , 
at the z'-th pointing. 
3. Using the value of Kq estimated in the previous point, we 



4 K 

extract a new set of gains K { ' with the equation 



Kf K = K x 



v: 



ef \ 



V ref ,0, 



(8) 



The procedure can be modelled by the following GNU R 3 
code: 

data<-data. frame (gain = iterative. dipole. gains, 

vref = 2 - signal. 4K/mean(signal. 4K) , 
dipole = dipole. amplitude. KCMB) 
fit<-lm(gain " vref + 8, data, weights = dipole) 
gains. 4K <- fit$coefficients[l] * data$dvref 

where iterative . dipole . gains, signal. 4K and 
dipole. amplitude. KCMB are three vectors containing 
the iterative gains Kf er before the smoothing filter, the 4K 
reference load signal VJ ef averaged over each pointing period, 
and the values of ATj (eq. 5). 

Unlike the procedure in section 7.1, in this case there is no 
need to smooth the stream of gains, as they share the same sta- 
bility of the voltages V[ ef . 

7.3. Colour Correction 

Colour corrections at intermediate spectral indices may be de- 
rived accurately from a quadratic fit to the values in Table 8. 
In addition, the data release includes the UcCC IDL package, in 
common between the two instruments, (Planck Collaboration IX 
2013) which calculates colour corrections and unit conversions 
using the band-averaged bandpass stored in the reduced instru- 
ment model (RIMO) file also included in the data release. The 
way this table has been computed is fully described in Planck 
Collaboration V (2013). 



8. Noise estimation 

The estimation of noise properties is fundamental in several as- 
pects of the Planck LFI data analysis. For instance, such mea- 
surements are used in the Monte Carlo simulations of noise nec- 
essary for power spectra estimations, as well as to determine 



3 http://www.r-project.org/. 
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Table 8. Multiplicative colour corrections cc(ff) for individual LFI Radiometer Chain Assemblies and for the band average maps. 



Horn Spectral index a 
-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

LFI-18 0.948 0.961 0.972 0.981 0.988 0.994 0.997 0.998 0.997 0.995 0.990 0.983 0.975 

LFI-19 0.856 0.878 0.899 0.919 0.939 0.957 0.975 0.991 1.006 1.020 1.032 1.043 1.053 

LFI-20 0.889 0.908 0.925 0.941 0.956 0.970 0.983 0.994 1.003 1.011 1.018 1.023 1.027 

LFI-21 0.917 0.933 0.947 0.960 0.971 0.981 0.989 0.996 1.001 1.004 1.006 1.006 1.004 

LFI-22 1.024 1.026 1.027 1.026 1.023 1.018 1.011 1.003 0.993 0.982 0.969 0.955 0.940 

LFI-23 0.985 0.991 0.996 0.999 1.001 1.002 1.002 1.000 0.997 0.993 0.988 0.982 0.975 

70 GHz 0.938 0.951 0.963 0.973 0.982 0.988 0.994 0.997 0.999 0.999 0.998 0.995 0.991 

LFI-24 0.978 0.984 0.988 0.993 0.996 0.998 0.999 1.000 0.999 0.998 0.996 0.993 0.989 

LFI-25 0.967 0.974 0.980 0.985 0.990 0.994 0.996 0.999 1.000 1.000 1.000 0.999 0.997 

LFI-26 0.957 0.966 0.973 0.980 0.985 0.990 0.995 0.998 1.000 1.001 1.002 1.002 1.000 

44 GHz 0.968 0.975 0.981 0.986 0.990 0.994 0.997 0.999 1.000 1.000 0.999 0.998 0.995 

LFI-27 0.948 0.959 0.969 0.978 0.985 0.991 0.995 0.998 1.000 1.000 0.998 0.995 0.991 

LFI-28 0.946 0.958 0.968 0.977 0.985 0.991 0.996 0.998 1.000 0.999 0.997 0.993 0.988 

30 GHz 0.947 0.959 0.969 0.977 0.985 0.991 0.995 0.998 1.000 1.000 0.998 0.994 0.989 



proper horn weights to be employed during the map-making 
process. In addition, inspection of noise properties during mis- 
sion life time is of paramount importance in tracking possible 
variations and anomalies in instrument behaviour/performance. 
In this respect we have improved our noise estimation pipeline 
compared to that used in Zacchei et al. (2011), implementing a 
Monte Carlo Markov Chain approach for the extraction of basic 
noise parameters such as the white noise level, knee-frequency 
and slope of the l/f noise at low frequencies. This allows for 
an un-biased estimation of the last two noise parameters com- 
pared to the previously implemented log-periodogram approach 
(Zacchei et al. 2011). 

The code works as follows. It assumes a functional form for 
noise spectrum: 



P(f) = o 2 



(9) 



with three basic parameters. It is also possible to work with a 
functional form with two more parameters: 

P(f) = cr 2 1 + + [J- 



(10) 



The latter could be useful when there are clearly two different 
behaviour in the low-frequency part of the spectrum where, be- 
side usual radiometric l/f noise, the signature of noise induced 
by thermal fluctuations appears. 

As for the white noise part, this is, as before, computed 
making a simple average of noise spectrum on the last 10% of 
frequency bins (the high-frequency tail of the noise spectrum). 
This percentage works well for almost all radiometers at 44 and 
70 GHz but is quite delicate for the 30 GHz radiometers which 
show typical values of knee-frequency around lOOmHz, and 
therefore require a smaller percentage to get an un-biased white 
noise estimation. Once white noise is computed, the code creates 
Markov Chains for the other parameters. Discarding the burn-in 
period of the chains we can directly get the expected value and 
variance of each noise parameter sampled from the chain distri- 
bution. 

The upper panel of Fig. 4 shows a typical spectrum at 
70 GHz with the simple log-periodogram fit (red line) and the 
new MCMC derived spectrum (blue line) superimposed. The 
right panel instead shows the distribution for knee-frequency and 
slope derived from the example spectrum. 



8. 1 . Updated noise properties 



Radiometer noise properties have been evaluated with the new 
code implementing the Monte Carlo Markov Chain approach 
described in the section above. As for the previous release of 
LFI data, we select radiometer calibrated data in chunks of five 
ODs and compute noise spectra with the roma iterative gener- 
alized least square (IGLS) map-making algorithm (Natoli et al. 
2001 ; de Gasperis et al. 2005), which includes a noise estimation 
tool based on the iterative approach described in Prunet et al. 
(2001). The actual algorithm is described therein and the inter- 
ested reader can find its application to LFI data in Zacchei et al. 
(2011). The output is a frequency spectrum on which the new 
MCMC code is applied. Results at radiometer level on white 
noise sensitivity are reported in Table 9, while Table 10 shows 
1 // noise parameters. These are computed taking the median of 
the five estimations made for different ranges of OD over the 
nominal mission time. 

As a matter of fact, time variations of the noise properties 
provide a valuable diagnostic of possible changes of instrument 
behavior. We know that there was quite a dramatic change in 
LFI operations, namely the switch over between the two sorp- 
tion coolers. We expect to see variations in the noise properties 
during the degradation process of the first cooler and at the onset 
of the redundant one. 

Figure 5 shows a set of noise frequency spectra for three LFI 
radiometers (LFI28M, LFI24S and LFI18M) based on the pe- 
riod considered for the 2013 data release. Some comments are 
in order. First of all, the white noise level is extremely stable 
in all three cases shown; this is also true for all the LFI ra- 
diometers. Also the knee-frequency and slopes are quite stable, 
as clearly demonstrated by the noise spectra until OD 326. After 
that time the spectra show an increasing noise and two slopes 
for the low-frequency part. The latter becomes more evident for 
spectra around OD 366 and OD 466, when the first cooler starts 
to be less effective and produces low-frequency thermal noise. 
This behavior is present at some level in all radiometers, but with 
different trends, ranging from the small effect shown by LFI24S 
to more prominent effects as shown by LFI28M and LFI18M. 
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Fig. 4. Left : typical noise spectrum at 70 GHz with simple log-periodogram fit (red line) and the new MCMC fit (blue line). Centre: 
sample distrubution for the knee-frequency and slope (Right) for the example spectrum. 
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Fig. 5. Time behavior of noise spectra for three LFI radiometers (Left LFI 18M 70 GHz, Centre LFI 24S 44 GHz and Right LFI 
28m 30 GHz). In all cases white noise is extremely stable as is the 1 // noise until OD 326. After that period, the degradation of the 
sorption cooler and the switch-over to the redundant cooler introduce higher thermal noise at the lowest frequencies. 



Table 9. White noise sensitivities for the LFI radiometers. 



White Noise Sensitivity 
Radiometer M Radiometer S 

[yuK CMB S 1 ' 2 ] [^K CMB S 1/2 ] 

70 GHz 



LFI18 51 1.7 ±1.7 466.3+1.7 

LFI19 579.8 ±1.6 554.1 ±1.6 

LFI20 587.5 ±1.4 619.7 ±2.1 

LFI21 451.6 ±1.7 560.9 ±1.7 

LFI22 489.9±1.5 531.0+2.1 

LFI23 503.4 ±1.8 538.8 ±1.9 

44 GHz 

LFI24 461.0 ±1.3 398.2 ±1.3 

LFI25 413.5 ±1.5 393.3 ±3.0 

LFI26 480.8 ±1.5 419.1 ±1.9 

30 GHz 

LFI27 282.2 ±2.1 304.7 ±2.0 

LFI28 318.2 ±1.9 286.8 ±2.1 



9. Map-Making 

The mapmaking pipeline was already described in detail in 
the Zacchei et al. (2011) and other papers. Here we will 
only describe the overall process reporting significant updates. 
Subsections will be dedicated to special processes such as the 
calculation of the noise covariance matrices, half-ring maps, 



survey-maps and Monte Carlo maps simulations used in the val- 
idation process. 

9.1. MADAM pipeline for frequency maps 

LFI frequency maps were produced by the Madam mapmak- 
ing code (Keihanen et al. 2010). The code takes as input the 
calibrated TOD for each radiometer, and the corresponding 
radiometer pointing data, in the form of three Euler angles 
(9, (f>, if/). The primary output consists of three pixelized Stokes 
maps (T,Q,U), which represent the temperature and polarization 
anisotropics of the observed sky. 

The algorithm is based on the destriping technique, where 
the correlated noise component is modelled by a sequence of 
offsets, or baselines. The amplitudes of these baselines are deter- 
mined through maximum-likelihood analysis. Higher-frequency 
noise, which is not captured by the baseline model, is assumed 
to be purely white. 

The noise model can be written as 

n — Fa + n (11) 

Here n' is the total noise stream, n represents white noise, a is 
a vector consisting of the baselines, and F is a matrix, consist- 
ing of ones and zeros, which spreads the baselines into a time- 
ordered data stream. 

Unlike conventional destriping, Madam also uses information 
on the known noise properties, in the form of a noise prior. This 
allows one to extend the destriping approach to very short base- 
line lengths, well below the scanning period of a minute. 
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Table 10. Knee frequency and slope (power of frequency) for the LFI radiometers. 



Knee Frequency 



Slope 



Radiometer M [mHz] 



Radiometer S [mHz] 



Radiometer M 



Radiometer S 



70 GHz 

LFI 18 
LFI 19 
LFI20 
LFI21 
LFI22 
LFI23 

44 GHz 

LFI24 
LFI25 
LFI26 

30 GHz 

LFI27 
LFI28 



15.3 +2.8 
11.9 ±1.3 
8.4 ±1.9 
39.3 ±4.0 

10.1 ±2.1 

30.2 ±1.4 



26.9 ±1.2 
20.1 ±0.6 
64.4 ±2.0 



175.1 : 
127.9 : 



:2.2 
:3.8 



18.3 ±1.6 
14.6 ±1.1 
6.0 ±1.7 
14.0 ±2.4 
15.9 ±7.1 
58.8 ±9.0 



73.0 ±7.9 

46.1 ±1.8 
43.8 ±8.9 



109.6 : 

43.9: 



:2.3 
:2.2 



-1.07 ±0.11 
-1.22 ±0.30 
-1.31 ±0.40 
-1.26 ±0.09 
-1.53 ±0.34 
-1.07 ±0.03 



-0.94 ±0.01 
-0.85 ±0.01 
-0.92 ±0.01 



-0.93 : 
-0.93 : 



:0.01 
:0.01 



■1.20 ±0.15 
-1.12 ±0.16 
-1.34 ±0.47 
-1.24 ±0.11 
-1.20 ±0.36 
-1.21 ±0.05 



-0.91 ±0.01 
-0.90 ±0.01 
-0.88 ±0.06 



-0.91 
-0.91 



:0.01 
:0.02 



The baseline length is a key parameter in the destriping 
technique. We chose the baseline length so that it corresponds 
roughly to one second of sampling time, while also requiring 
it to be an integral number of samples. The chosen baseline 
lengths were 33, 47, and 79 samples, for the 30, 44, and 70 GHz 
channels, respectively. These correspond to 1.0151, 1.0098, and 
1 .0029 s, respectively. 

The baseline length of one second is a reasonable compro- 
mise between computational burden and the quality of the final 
map. Shortening the baseline further, below one second, has very 
little effect on the residual noise. 

When mapmaking is run with a short baseline, the noise prior 
plays an important role. There are not enough crossing points be- 
tween the one-second data sections to determine the baselines 
without inserting additional a priori information to constrain 
the solution. In this case, the required additional information is 
the noise prior, which was constructed from the parametrized 
noise model (knee frequency, slope, white noise variance), based 
on the parameters given in Tables 9 and 10. From these we 
compute the expected covariance between the noise baselines, 
C a = (flfl T ). The exact derivation is given in Keihanen et al. 
(2010). Another important role of the noise prior is to suppress 
the signal error, which would otherwise increase rapidly with 
decreasing baseline length. Flagged data sections were handled 
by setting the white noise variance formally to infinity for those 
samples, but not altering the baseline pattern. The flagged sam- 
ples thus do not contribute to the final map in any other way 
than serving as space-holders which keep the correct distance 
between unflagged samples on both sides of the gap. This is es- 
sential for the noise prior to be applied correctly. 

The use of a priori information has the danger of hiding 
problems in the data, since the prior may drive the solution to 
a correct-looking result, even if the data alone were not in com- 
plete agreement with it. To avoid this pitfall, we computed for 
comparison a subset of the maps with a one-minute baseline, 
without using a noise prior, and compared the maps visually. No 
artifacts could be seen in either of the maps. Despite the trivial- 
ity of the check, similar sanity checks helped to reveal bugs at 
earlier stages of the development of the data processing pipeline. 

Madam uses HEALPix pixelisation to discretize the sky. Each 
data sample is assigned as a whole into the pixel where the cen- 
ter of the beam falls. We used resolution = 1024 for all LFI 



frequencies. The average width of one pixel at this resolution is 
3.5'. The chosen resolution gives complete sky coverage for the 
70 GHz channel. At 30 GHz and 44 GHz, where the sampling 
frequency is lower, there remain individual unobserved pixels, 
which are marked by a special value in the product maps. In 
the 30 GHz frequency map there are 158 and and in the 44 GHz 
frequency map 250 missing pixels (0.0013% and 0.0020%, re- 
spectively). These count also the pixels which were covered with 
insufficient polarisation angle coverage, and for which we could 
not recover the polarisation component. The pixel size is to be 
compared with the beam resolution at each frequency (Table 3). 

The maximum-likelihood analysis, which lies behind our 
mapmaking algorithm, and derivation of the destriping solution, 
are presented in Keihanen et al. (2010). Here we quote the most 
important formulas for easier reference. 

The vector of baseline a is solved from the linear equation 

(F T C- l ZF + C- 1 )a = F r C- 1 Zy, (12) 

where 

z = i-(p t c- [ p) 1 p t c- 1 . (13) 

The final map is then constructed as 

m = (P T C; i 'Py'P T C-\y - Fa), (14) 

where a is the baseline solution. Here F is a matrix, consisting 
of ones and zeros, which spreads the baselines into time-ordered 
data. P is the pointing matrix, which picks values from the T,Q,U 
maps and spreads them into time-ordered data, and y is the ob- 
served data stream. The covariance matrices C„ and C w represent 
the a priori known properties of the baselines and the white noise 
component, respectively. The white noise component is assumed 
to be uncorrelated, but not necessarily with uniform variance. 
Matrix C w is thus diagonal, but not constant. 

Matrices P, F, C w are very large, but sparse. They are not 
constructed explicitly, but the operations represented formally as 
matrix multiplication above, are performed algorithmically. For 
instance, multiplication by P T represents an operation where the 
time-ordered data are coadded on a sky map. 

We solve Eq. 12 through conjugate gradient iteration. 
Convergence was reached typically after 20-100 iterations, de- 
pending on the sky coverage and radiometer combination. Full 
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mission maps, where the whole sky was covered, typically con- 
verge faster than single-survey maps. 
The observed signal can be written as 

yi = T(fi>d + Q{ud cos(2^) + U{iod sin(2^), (15) 

where w, is the sky pixel where the sample is assigned to, and 
if/i defines the beam orientation. The elements of the pointing 
matrix P consist of ones and cosine and sine factors picked from 
this equation. 

When constructing single-horn maps, we included only the 
temperature component of the sky into the computation. In the 
case of frequency or horn-pair maps we included three map 
components, corresponding to the I,Q, U Stokes components, al- 
though only the I component maps are included in the 2013 data 
release. 

We deviate from the formulation of the original Madam paper 
in that we have written C w in place of C„. This reflects the fact 
that the C K matrix we insert into the solution is not necessar- 
ily the white noise covariance, but rather a user-defined weight- 
ing factor. We decided to diverge from the exact maximum- 
likelihood solution in order to have better control over system- 
atics. The weight was constructed according to a horn-uniform 
weighting scheme, as follows. The weight for a given horn was 
chosen to be 



where <tm and cr? are the white noise sensitivities of the two 
radiometers of the horn, computed from the values given in Table 
3. The weights were identical for radiometers of same horn. The 
formula above is applied for non-flagged samples. For flagged 
samples we set C" 1 = 0. 

When horn-uniform weighting is applied, the polarization 
maps become dependent solely on the signal difference between 
M and S radiometers, apart from a small leakage due to the 
fact the polarization sensitivities are not exactly at 90°from each 
other. Many systematic effects, which are equal or strongly cor- 
related within a horn, cancel out in the differentiation. In par- 
ticular, the horn-uniform weighting has the benefit of reducing 
spurious polarization signals arising from beam shape mismatch, 
since the beam shapes of a radiometer pair sharing a hom are 
typically quite similar (though not identical). Complete cancel- 
lation also requires that the same flags are applied to both data 
streams. Therefore, if a sample for one radiometer was flagged, 
we discarded the corresponding sample for the other radiometer 
as well. 

The covariance of residual white noise in the map solution is 
obtained from 

c wn = {P T c- l py 1 p T c- 1 c n c- 1 P{P T c- l P)- x (17) 

where now C„ is the actual white noise covariance of the time- 
ordered data. For a given C„, the residual noise is minimized 
when C w = C„, in which case C wn = (P T C~ 1 P)~ l . In using a 
different weighting we thus sacrifice some noise sensitivity for 
better removal of systematics. 

Note that altering matrix C w does not introduce bias into the 
solution, as may easily be verified by inserting y = Pm in the 
solution above. It merely affects the level of residual noise. 

The white noise covariance only takes into account the un- 
correlated component of the noise. The computation of a full 
noise covariance matrix, which also captures the residual corre- 
lated noise, is discussed in a separate section below. 

The destriping solution assumes that the sky signal is uni- 
form within one pixel. This is not strictly true, which gives rise 



to the signal error: signal differences within a pixel are falsely 
interpreted as noise, which leads to spurious striping, especially 
in the vicinity of point sources or in other regions where the sig- 
nal gradient is large. Most of the effect comes from the Galactic 
region. 

Similarly, mismatch in frequency response between ra- 
diometers gives rise to spurious striping, as different radiome- 
ters record slightly different signals from the same source. In 
this case also the main effect comes from the Galactic region, 
where foreground emission is strong. 

To reduce these undesired effects, we masked out the 
Galactic region and point sources in the destriping process. We 
used crossing points of rings only outside the masked region to 
solve the noise baselines, though all the data was included in the 
construction of the final map product. The masks for 30, 44, and 
70 GHz left uncovered 78.7%, 89.4%, and 89.7% of the sky, re- 
spectively (fraction of sky included in the analysis). Some of the 
baselines fall completely inside the masked region, but can still 
be recovered with reasonable accuracy with help of combined 
information of the neighbouring baselines and the noise prior. 

In Table 1 1 we list the delivered maps. All the maps have 
HEALPix resolution N side = 1024. 

9.2. Low Resolution data set 

To fully exploit the information contained in the large scale 
structure of the microwave sky, pixel-pixel covariances are 
needed in the maximum likelihood estimation of the CMB power 
spectrum. However, full covariance matrices are impossible to 
employ at the native map resolution, because of resource lim- 
itations. A low-resolution dataset is therefore required for the 
low-^ analysis. This dataset consists of low -resolution maps, and 
descriptions of residual noise present in those maps given by 
pixel-pixel noise covariance matrices (NCVMs). 

The low-resolution dataset can currently be utilized effi- 
ciently only at resolution N s id e = 16, or lower. All the low- 
resolution data products are produced at this target resolution. 

We will first discuss the production of the low-resolution 
maps, and then continue the discussion on the NCVMs. 

9.2.1. Low-resolution maps 

Different schemes to produce the low-resolution maps are dis- 
cussed in Keskitalo et al. (2010). We chose to construct low- 
resolution maps by downgrading the high-resolution maps to the 
target resolution, N s id e = 16, using a noise- weighted downgrad- 
ing scheme. A low-resolution map was obtained by applying to 
its high-resolution counterpart the operation 

mi = {P]C^PxY l X(,PlC- l P b )m h = Dm h , 

where X qp sums the high-resolution pixels to the low-resolution, 

{1 , p subpixel of q 
0, otherwise. 

Here subscripts h, 1 refer to the high {N^ = 1024) and low 
(A^side = 16) resolution versions of the pointing matrix. The same 
matrix X downgrades the pointing matrix, P\ = PhX T . The re- 
sulting map is identical to the one we would get by solving the 
baselines at the higher resolution, but then binning the map to a 
lower resolution. 

By using the noise-weighted downgrading scheme, we get 
an adequate control over signal and noise in the resulting map. 
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Map Horns OD range Baseline (s) Sky coverage (%) 



30 GHz nominal 27,28 91-563 1.0151 99.999 

30 GHz survey 1 27,28 91-270 1.0151 97.205 

30 GHz survey 2 27,28 270-456 1.0151 97.484 

44 GHz nominal 24,25,26 91-563 1.0098 99.998 

44 GHz survey 1 24, 25, 26 91-270 1.0098 93.934 

44 GHz survey 2 24,25,26 270-456 1.0098 93.310 

70 GHz nominal 18,19,20,21,22,23 91-563 1.0029 100.000 

70 GHz survey 1 18, 19, 20, 21, 22, 23 91-270 1.0029 97.938 

70 GHz survey 2 18, 19, 20, 21, 22, 23 270-456 1.0029 97.474 



If we were to calculate a low-resolution map directly at the tar- 
get resolution, the signal error would be larger due to sub pixel 
structures, while the noise level would be lower. 

After downgrading, the temperature component was 
smoothed with a symmetric Gaussian window function with 
FWHM = 440', while the polarization components were left 
unsmoothed. Smoothing was applied to alleviate aliasing due to 
high frequency power in the map. The polarization components 
need to be treated differently, since the cosmological informa- 
tion contained in polarization has several orders of magnitude 
poorer signal-to-noise ratio. Compared to the approach proposed 
in Keskitalo et al. (2010), we intentionally changed the order 
of downgrading and smoothing to better deal with noise. The 
aliased power is negligible at scales unaffected by the smooth- 
ing operator. 

9.2.2. Noise covariance matrices 

The statistical description of the residual noise present in a low- 
resolution map is given in the form of a pixel-pixel noise co- 
variance matrix, as described in Keskitalo et al. (2010). We must 
apply the processing steps chosen in the map downgrading to 
NCVMs for consistency. Some approximations were however 
required: some are inherent to the Madam map-maker NCVM 
module, while some were made due to performance issues. 

The pixel-pixel noise covariance matrix for generalized de- 
striping is 

N = (P T (C w +FC a F T r 1 P)-\ 
which can be written in a dimensionally reduced form as 

AT 1 = P T C- l P - P T C- l F(F T C- l F + C- l y l F T C; i l P (18) 

Applying Eq. 1 8 in practice, lead to the need of inverting a sym- 
metric 3Afpj x x 3A^pi x matrix in a later analysis step. Because the 
inverse NCVMs are additive, we divided the computations into 
a number of small chunks to save computational resources. We 
first calculated, using Eq. 18 one inverse NCVM per radiometer 
per survey at the highest possible resolution N s id e = 32 permit- 
ted by computer resources. Later we combined the individual 
inverse matrices to the actual inverse matrices. 

To obtain the noise covariance from its inverse, the matri- 
ces were inverted using the eigendecomposition of a matrix. The 
monopole of the temperature map cannot be resolved by the 
map-maker, and thus the matrix becomes singular. Therefore the 
ill-determined mode was left out of the analysis. 

These intermediate resolution matrices are then downgraded 
to the target resolution. The downgrading operator is the same D 
as for the map downgrading, but Pj, is replaced with P„ i.e., with 
the intermediate resolution pointing matrix. The downgraded 



matrix is 

N, = DND T . (19) 

As a final step the same smoothing operator was applied to 
the temperature component of the matrices as was applied to the 
low-resolution maps. 

The noise covariance matrices were calculated with two dif- 
ferent sets of noise parameters. One set covered the entire mis- 
sion under consideration (values given in Tables 9 and 10), 
while the other had individual noise parameters for each survey. 

Equation 18 describes the noise correlations in a situation 
where the noise baselines are solved at the resolution of the fi- 
nal map. For an exact description of the correlations, we should 
construct the matrix at resolution A^ s ;d e = 1024, and downgrade 
the inverted matrix to the target resolution N s id e = 16. This is not 
feasible, due to the size of the matrix. We therefore construct the 
matrix at the highest possible resolution = 32, and down- 
grade it to the target resolution. 

The same formula inherently assumes that individual detec- 
tors are weighted according to their white noise levels, as sug- 
gested by maximum likelihood analysis. In map-making, how- 
ever, we apply horn-uniform weighting, to have better control 
over systematics, as explained in the previous section. Also, the 
formulation does not take into account the effect of the destrip- 
ing mask. As results of these idealizations, the covariance matrix 
is an approximate description of the residual noise correlations. 
We have performed^ 2 tests to assess the effect of each idealiza- 
tion individually. 

Figure 6 illustrates the effect of removing the approximations 
inherent in the NCVM computation: no horn-uniform weighting; 
the destriping resolution, Nueaa., equalling the map resolution; 
and no masking in the destriping phase. Each non-ideal factor 
increases the discrepancy seen in the x 1 test. 

We chose to use shorter baselines in the NCVM production 
than in the map-making. They were 0.25 s / 8 samples, 0.5 s / 24 
samples, and 0.5 s / 39 samples for 30 GHz, 44 GHz and 70 GHz, 
respectively. Since the knee frequencies were higher than an- 
ticipated prior to launch, short baselines model the noise better 
(Keskitalo et al. 2010). We additionally found out that reducing 
baseline length in the NCVM calculation affects the x 2 statis- 
tics more than changing baseline length in the map-making (see 
Fig. 7). 

9.3. Half ring noise maps 

In order to estimate the noise directly at the map level and in 
the angular power spectra, we produced half-ring maps (hi and 
h.2) with the same pipeline as described in Sect. 9.1, but using 
data only from the first or the second half of each stable point- 
ing period. These half-ring maps contain the same sky signal, 
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=1024, Destr. mask 



N D .„ =1 024, Destr. mask, 
Horn-uniform weighting 



Chi" 

Fig. 6. Reduced^- 2 statistics from 25 noise only maps for 30 GHz 
2013 delivery. The NCVM was calculated using 0.25 s base- 
lines, while the simulations were made with 0.5 s baselines. The 
number of idealizations in the noise only simulations decreases 
from top to bottom. The first set of simulations (plotted in black) 
contains the same approximations that are made in the NCVM 
calculation. The last set of simulations (plotted in purple) cor- 
responds to the standard map-making options: the horns are 
weighted uniformly; destriping resolution, Abestr., is 1024; and 
a destriping mask is applied. 



since they result from the same scanning pattern on the sky. 
Therefore the difference of maps hi and I12 captures any noise 
that varies faster than half of the duration of the pointing period, 
i.e., the noise whose frequency is / > 1/20 min = 0.85 mHz. 
The procedure of calculating the half-ring maps and their hit 
count weighted difference maps is described in more detail by 
Zacchei et al. (201 1) and Planck Collaboration ES (2013). 

The use of half -ring maps in the validation of data and noise 
estimates will be explained in Sect. 12. In addition, the half-ring 
maps were an integral part of the component separation process 
Planck Collaboration XII (2013) and likelihood codes Planck 
Collaboration XV (2013). 

9.4. Noise Monte Carlo simulations 

Simulated noise timelines were produced according to the three- 
parameter noise model (white noise sensitivity, knee frequency, 
slope), using the estimated values for these parameters given in 
Tables 9 and 10. Maps were made from this simulated noise 
using reconstructed flight pointing and the same Madam param- 
eter settings as were used for the flight maps. These steps were 
repeated to produce a Monte Carlo set of 1000 realizations of 
noise maps for different radiometer combinations, including fre- 
quency maps and 70 GHz horn-pair maps. 

For LFI, this Monte Carlo (MC) work was done in two 
stages, with two partially different pipelines, first an LFI MC 
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Fig. 7. The reduced x 2 statistics from noise only maps for 
30 GHz 2013 delivery. Upper, the noise only simulation set is 
fixed, while the NCVM baseline length changes. The baseline 
lengths in the NCVM calculation were 1 s, 0.5 s, and 0.25 s, 
and are plotted in black, blue, and red respectively. Lower, the 
NCVM is fixed, while the noise only simulation varies. Two 
baseline lengths were chosen 1 s (plotted in black) and 0.25 s 
(red). 



in close connection with the map-making from the LFI flight 
data, and then as a part of the joint LFI/HFI full-focal plane 
"FFP6" simulations. This work was divided between two su- 
percomputing centers, CSC-IT Center for Science in Finland 
and the National Energy Research Scientific Computing Center 
(NERSC) in U.S. 

The noise MC maps provide a statistical distribution of noise 
maps that can be compared to the half -ring noise maps (Sec. 9.3), 
to see how well maps from the noise model match the real flight 
noise in the half-ring noise maps. We note, however, that the 
half -ring noise maps do not represent properly the noise in the 
flight maps for the lowest frequencies, i.e., for time scales of the 
order of half the pointing period or longer. 

For low-resolution studies the maps were downgraded to 
Miide = 32 and N s n e = 16 HEALPix resolution, using the same 
procedure as for the flight maps. These can be compared to the 
low-resolution noise covariance matrices discussed in Sec. 9.2, 
which were generated from the same noise model, but are based 
on some approximations. Thus this comparison, see Fig. 6 and 
Fig. 7, shows the effect of these approximations on the NCVM. 

In addition the noise MC maps were used in power spectrum 
estimation, component separation (Planck Collaboration XII 
2013), and in non-Gaussianity estimation (Planck Collaboration 
XXIII 2013) and (Planck Collaboration XXIII 2013). 

9.5. Overview of LFI map properties 

Figures 8 to 10 show the sky frequency maps created from LFI 
data. The top row of each figure shows the temperature maps 
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based on the nominal mission data, while the second row shows 
the differences between maps, « m , made of the first and second 
half of each stable pointing period (half-ring maps) weighted by 
the hit count calculated from equation 21 . These maps provide 
an intuitive evidence of the level and distribution of the residuals 
and are calculated as 



hi - h 2 

Whit 



where the hit count weight is 



Whit = a hit 



full 



1 1 

hit; hit2 



(20) 



(21) 



Here hitf u n = hitj + hit2 is the hit count of the full map m, while 
hiti and hit2 are the hit counts of the haf-ring maps hi and /12, 
respectively. 

The third row shows the difference in the sky between the 
first and second survey, and gives information on longer time 
scale variations. In the half-ring difference maps a darker stripe, 
corresponding to two observation days, is visible due to higher 
noise (it is barely visible also in the frequency maps). This is due 
to the fact that in the first days of observation the LFI instrument 
was affected by an occasional bit-flip change in the gain-setting 
circuit of the data acquisition electronics, probably due to cosmic 
ray hits; at that time we did not know the cause of the problem 
and we were obliged to perform some operations to identify and 
solve it. The data acquired during this period, and contributing 
to the visible ring, have been flagged out. One other clear feature 
can be seen on the Galactic plane in the survey difference maps, 
especially at 30 GHz. This apparent split in intensity is due to 
the beam ellipticity: the elliptical beam had a different orienta- 
tion relative to the Galaxy in the first survey compared with the 
second. 



10. Polarization 

The LFI data processing has included analysis of polarization 
from the beginning, but polarization results are not included in 
the current data release and scientific analysis, because the level 
of systematic errors in the maps remains above acceptable levels. 
In this section we briefly outline the polarization-specific steps 
in the data analysis, quantify the residual systematics, and sketch 
how we expect to correct them for the next data release. 

To an excellent approximation (see Sect. 9.1) the LFI po- 
larization maps (Stokes parameters Q and U) are derived from 
the difference between the calibrated signals from the two ra- 
diometers in each RCA, so-called main-arm (M) and side-arm 
(S), which are sensitive to orthogonal polarizations. Any differ- 
ential calibration errors between M and S cause leakage of total 
intensity into the polarization maps. Such mismatch arises from 
three main causes: 

- differences between the beam profiles of M and S; 

- errors in the gain calibration; 

- differential colour corrections between M and S due to dif- 
ferences in their bandpasses (the 'bandpass leakage') effect. 

All these effects are described in Leahy et al. (2010). Only the 
bandpass leakage requires a special step in the calibration, al- 
though polarization imposes stringent requirements on the accu- 
racy of gain calibration which have driven our choice of calibra- 
tion scheme (Planck Collaboration V 2013). The control of sys- 
tematics for polarization also require precise cancellation of the 



M and S arm signals, which underlies our choice to use identical 
pointings for the data from the two radiometers in each hom, de- 
spite a small amount of beam squint between the polarizations, 
and also the decision to use horn-uniform weighting. Neither of 
these choices results in significant degradation to the total inten- 
sity maps although they are slightly non-optimal. 

The principle instrumental factor controlling bandpass leak- 
age is the effective frequency mismatch between M and S de- 
tectors, a = (vs - vm)/2vo, which must be combined with es- 
timates of the "leakage amplitude" of the foreground emission, 
L = (jSfg - /?cmb)7V Because the leakage amplitude relies on 
products from component separation, maps of L are available at 
no higher resolution than that of the 30GHz channel (33.16 '), 
and are most reliable at lower resolution ( 1 °), where our analysis 
can incorporate the WMAP 22 GHz maps. 

In principle the a-factors can be estimated from the bandpass 
profiles measured in the ground calibration campaign (Zonca 
et al. 2009), but as anticipated by Leahy et al. (2010), more ac- 
curate values can be found from the flight data; a detailed de- 
scription of our approach to this will be discussed in a future 
paper. We estimate that our a-factors are currently accurate to 
about 0.05%, based on the scatter in values derived from differ- 
ent calibrators. A bigger problem at present is accurate evalu- 
ation of the leakage amplitude, which requires not only excel- 
lent separation of CMB and foreground emission but also accu- 
rate estimates of the foreground spectral index within the band. 
Currently, the combined uncertainty in the bandpass leakage cor- 
rection is about 0.3% of the local foreground intensity. This is 
comparable to the mean polarization fraction along the Galactic 
plane. 

One other parameter in principle should be calibrated for po - 
larization: the precise orientation of the polarization response for 
each feed horn. As noted by Leahy et al. (2010), knowledge of 
horn orientation was expected to be good to better than 1 ° from 
the construction tolerances, even though it was not possible to 
calibrate them directly. In flight, we check these values by ob- 
servations of the Crab nebula, and the results confirm the orien- 
tations to within a few degrees, which is sufficient for analysis of 
the £-mode spectra. At a higher level of precision our estimates 
of the response orientations are still affected by the uncertainty 
in a-factors, and so as yet we have no evidence to reject the nom- 
inal orientation angles. 

The most interesting cosmological signal visible in LFI po- 
larization is the very large-scale (£ < 10) E'-mode peak due 
to reionisation, at a typical brightness level of 0.3 fiK. Crucial 
tests of the reliability of this signal are that the B-mode and 
EB cross-correlation spectra should contain negligible signal, 
since a cosmological fi-mode signal at this level would corre- 
spond to a tensor-to-scalar ratio significantly larger than cur- 
rent upper limits, while the cosmological EB mode is precisely 
zero. The primary channel for cosmology is 70 GHz, which has 
the least foreground contamination of all Planck channels; our 
likelihood pipeline estimates these spectra using a conservative 
Galactic mask and corrects for residual foregrounds based on 
Planck 30 GHz map and/or WMAP maps. With our current cal- 
ibration, both these spectra contain residuals at a level compa- 
rable to the expected E-mode signal from reionisation. Hence, 
although the latter is apparently detected, we cannot be confi- 
dent that the signal is real. This situation is better illustrated in 
Fig. 12 where we report, at all three LFI frequencies, null tests 
spectra from survey-survey differences. A comment is in order 
here. At all frequencies the null test EE spectra are in very good 
agreement with the expected noise level as traced by half-ring 
difference maps at multipoles larger of few tens: this is an indi- 
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Fig. 8. LFI maps at 30 GHz. The first row gives the intensity expressed in (iKcmb- The second row shows the difference between 
maps made of the first and the second half of each stable pointing period. The third row display the difference between the first and 
second survey maps. 
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Fig. 9. LFI maps at 44 GHz. The first row gives the intensity expressed in ^Kcmb- The second row shows the difference between 
maps made of the first and the second half of each stable pointing period. The third row display the difference between the first and 
second survey maps. 
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Fig. 10. LFI maps at 70 GHz. The first row gives the intensity expressed in ^uKcmb- The second row shows the difference between 
maps made of the first and the second half of each stable pointing period. The third row display the difference between the first and 
second survey maps. 
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cation of the data quality at these multipoles. However at very 
low-^ residuals are present especially at 30 GHz and at 70 GHz 
and these unable us for a proper characterisation of the cosmo- 
logical signal. 

We have simulated the impact of numerous systematic er- 
rors to see if they can explain the observed residuals, including 
foreground correction, bandpass mismatch, Galactic straylight 
(i.e., leakage through the far sidelobes), and gain errors. None of 
these simulations have individually generated artefacts as large 
as those observed. The most likely candidate seems to be in the 
combination of far sidelobe correction and calibration errors. As 
described in Planck Collaboration V (2013), uncertainty in the 
far sidelobe pattern is one of the dominant contributors to our 
calibration uncertainty, as well as making our estimates of the 
additive effect of Galactic straylight quite uncertain. 

1 1 . Power Spectra 

LFI temperature power spectra are computed from frequency 
maps using a cROMAster, an implementation of the pseudo-Cf 
method described in Hivon et al. (2002). We extend it to derive 
both auto- and cross-power spectra (see Polenta et al. (2005) 
for a comparison between the two estimators). Noise bias and 
covariance matrices have been computed through the full focal 
plane simulations version 6 "FFP6" that include 1000 realiza- 
tion of both signal and noise maps consistent with Planck data. 
The angular response of the instrument is accounted for by using 
the beam window functions presented in Planck Collaboration 
IV (2013). Coupling kernels to correct for uncompleted sky 
coverage are computed as described in Annex. B of Planck 
Collaboration XV (2013). We have masked the Galactic Plane 
and point sources using masks described in Sec. 3 of Planck 
Collaboration XII (2013). In particular, we have used a 70% 
Galactic mask for 44 and 70 GHz (leaving 70% of the sky un- 
covered), while we have used a 60% Galactic mask for 30 GHz. 

We show in Fig. 11 the 30, 44, and 70 GHz temperature 
power spectra. These have been produced from frequency maps 
without performing component separation. Nevertheless, there is 
a clear agreement between the observed spectra and the Planck 
likelihood code best fit curve (Planck Collaboration XV 2013) 
when we add a simple foreground component to account for un- 
resolved point source residuals. 

12. Data validation 

In order to assess and verify the quality of LFI data produced 
within the data-analysis pipeline, a set of null tests was prepared 
and performed. The main goal is to validate LFI science data 
by highlighting possible instrumental systematic effects. These 
include effects which are either properly corrected or accounted 
for in the pipeline; or are related to known changes in the opera- 
tional conditions of the LFI instrument (e.g., the switch over of 
sorption cooler) or to intrinsic instrument properties coupled to 
the sky like stray-light from sidelobes. Such tests are also useful 
for a detailed analysis of the processing steps implemented in 
the pipeline. Finally, such null tests may discover unanticipated 
problems (e.g., related to different calibration approaches). 

Null tests were carried out on blocks of data on different time 
scales ranging from pointing periods to one year of observation, 
and at different unit levels (radiometer, horn, horn-pairs within a 
given frequency and at frequency level) for both total intensity 
and, when applicable, for polarization. Such approach is quite 
demanding in terms of all possible combinations and specific 
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Fig. 11. Temperature power spectra at 30, 44, and 70 GHz. 
Dashed lines correspond to the Planck Likelihood Code best- 
fit plus a foreground component to account for unresolved point 
sources. 



tools to handle this work-load have been prepared by a dedicated 
team that creates a parallel code, responsible for both the actual 
computation of the null-tests and the creation of the output maps 
and spectra, and the visualization code that creates a report (a 
HTML page) with all the output from the null tests. 

The kind of effects probed by null tests depend on the com- 
bination of data and time scale treated. For example, differences 
at horn level between odd and even number surveys clearly re- 
veal the impact of sidelobes since the sky covered is exactly the 
same but beam orientation is not. On the other hand, differences 
between horns for the entire data period may reveal calibration 
issues and/or changes in the operational or instrumental condi- 
tions. 

12.1. Null tests summary 

It is important also to set a pass-fail criterion for such null tests. 
In general a failure of a specific test reveals some issue in the 
data or in the data analysis. These have to be carefully studied 
and specific actions have to be taken in order to at least miti- 
gate the non ideal result revealed by the test. A simple figure 
of merit would be the actual level of noise in the data derived 
from half-ring differences maps. Any departure from the half- 
ring difference maps noise level would be an indication of some 
problems. For example (Planck Collaboration III 2013) null test 
power spectra are used to check the total level of systematic ef- 
fects in the data. In Fig. 12 we report results at frequency level of 
survey difference null tests for both TT and EE spectra compared 
to the noise level derived from half-ring differences maps. 

It is worth mentioning that from such results we have care- 
fully analyzed our calibration pipeline and particularly the treat- 
ment of sidelobes. This has been corrected, updated with the 
inclusion of both the intermediate beam as well as the in-band 
beam behavior based on simulations. This will be the final ap- 
proach for data calibration in the next release. 

Although it is clear from these results that a proper treatment 
of sidelobes is necessary for final refinements in calibration, it 
is important to note that the overall amplitude of such effects is 
well below the CMB signal in total intensity, leaving the analysis 
of LFI temperature maps totally unaffected. 
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Fig. 12. Null tests results: power spectra from survey difference maps (SS) compared to half-ring difference maps (jk) noise level. 
Some excess at low multipoles is clearly visible at 30 GHz (Left), where the main source has been identified as sidelobe contribution. 
At 44 (Centre) and 70 GHz (Right), those null tests provide an indication of less impact of the low-f contribution both in TT and in 
EE although residuals are still present. For multipoles larger than few tens, null tests EE spectra follow the expected level of noise 
as traced by half -ring difference. 



12.2. Half ring test 

The middle panels of Figs. 8, 9, and 10 show the noise cal- 
culated from the hit-count-weighted half-ring difference maps at 
the native N s \d e = 1024 resolution, as described by Zacchei et al. 
(2011) and Planck Collaboration ES (2013). As a first quality 
check of the maps (and as one of the tests of the whole data 
processing pipeline up to the maps) we tested both numerically 
and visually that the noise maps divided pixel-by-pixel by the 
square root of the white noise covariance maps (see again Planck 
Collaboration ES (2013)) were approximately Gaussian with 
rms near to unity; the results were 1.0211, 1.0089, and 1.0007 
for 30, 44, and 70 GHz, respectively. 

The half-ring difference maps n m are the most direct mea- 
sure of the noise in the actual maps. The other noise estimates 
(NCVM and noise Monte Carlo) rely on specific modeling of 
the noise in the TOD, and this modeling can be validated by 
comparing the results to the half-ring difference maps. For this 
purpose we calculated the temperature and polarization (E and 
B mode) auto-correlation and cross-correlation angular power 
spectra of the noise by anafast using the half-ring difference 
maps and compared to these the results from the white noise co- 
variance matrices (WNCM) calculated by both madam and the 
noise Monte Carlo simulations (Sect. 9.4). Figure 13 illustrates 
such a comparison. Further, we calculated the mean C( for the 
high-£ tails (t — 1 150 — 1800) of the noise angular power spec- 
tra and took the ratio to the WNCM estimate; see Figure 14. As 
expected, there is some residual 1 // noise even in the high-^ re- 
gion, i.e., the full noise MC leads to slightly higher noise predic- 
tion than the WNCM or binned white noise from the noise MC. 
The residual l/f noise is of the order of 2.5% at 30 GHz, 1.0% 
at 44 GHz, and 0.1% at 70 GHz. Between the noise MC and 
the direct noise calculation from the half-ring difference maps 
we find good consistency: the high-/' noise from noise MC is 
only 0.8% higher at 30 GHz, 0.2% higher at 44 GHz and 0.1% 
lower at 70 GHz than the result from the half -ring difference. It 
should be noted that the error bars of the noise MC do not in- 
clude at this stage the uncertainty of noise parameters indicated 
in Tables 9 and 10. More such comparison are reported in Planck 
Collaboration ES (2013). 

12.3. Intra frequency consistency check 

We have tested the consistency between 30, 44, and 70 GHz 
maps by comparing the power spectra in the multipole range 



around the first acoustic peak. In order to do so, we have re- 
moved the estimated contribution from unresolved point source 
from the spectra presented in Sec. 1 1 . We have then built the 
scatter plots for the three frequency pairs, i.e., 70 vs 30 GHz, 
70 vs 44 GHz, and 44 vs 30 GHz, and performed a linear fit ac- 
counting for errors on both axis. 

The results reported in Fig. 15 show that the three power 
spectra are consistent within the errors. Moreover, please note 
that current error budget does not account for foreground re- 
moval, calibration, and window function uncertainties. Hence, 
the resulting agreement between spectra at different frequen- 
cies can reasonably be considered even more significant. We 
also compared the flux densities of compact sources at the three 
LFI frequencies, derived from the PCCS, Planck catalogue of 
compact sources (Planck Collaboration XXVIII 2013) and find 
these in acceptable agreement. These tests are described briefly 
in Planck Collaboration XI (2013). 

12.4. 70 GHz internal consistency check 

We use the Hausman test (Polenta et al. 2005) to assess the con- 
sistency of auto and cross spectral estimates at 70 GHz. We de- 
fine the statistic: 

H e = (C t - C e ) I ^Var{Cc-C c ), (22) 

where Q and Q represent auto- and cross-spectra, respectively. 
In order to combine information from different multipoles into a 
single quantity, we define the following quantity: 

, [Lr] 

B L (r) = — Y u H i ,re\Q,l} (23) 
yL f,_2 

where [.] denotes integer part. The distribution of Bi(r) con- 
verges (in a functional sense) to a Brownian motion process, 
which can be studied through the statistics s\ = sup r B/,(r), 

S2 = sup r \BL(r)\ and 53 = JT B 2 L (r)dr. Using the "FFP6" sim- 
ulations we derive the empirical distribution for all the three test 
statistics and we then compare them with the results obtained 
from Planck data itself (see Fig. 16). The Hausman test shows 
no statistically significant inconsistencies between the two spec- 
tral estimates. 

As a further test, we have estimated the temperature power 
spectrum for each of three horn-pair map, and we have compared 
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Fig. 13. Comparison of the noise angular power spectra of intensity maps: red from half -ring difference maps, blue from white noise 
covariance maps produced by madam, and black 16%, 50%, and 84% quantiles from noise MC for each C( . (Note that in the noise 
MC case no errors were propagated from Tables 9 and 10; only the median values of the three noise parameters were used.) 
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Fig. 14. The ratio of mean noise angular power at high- f ({ = 
1150 — 1800) to the white noise estimate from white noise co- 
variance matrices: red from the hit-count-weighted half-ring dif- 
ference, black from 101 full noise Monte Carlo maps, and blue 
the binned white noise from the noise Monte Carlo. The "ref- 
erence" white noise levels (green) are 9.8 x 10~ 15 K 2 , 13.1 x 
10~ 15 K 2 , and 10.1 x 10~ 15 K 2 , for 30, 44, and 70GHz, respec- 
tively. (Note that in the noise MC no errors were propagated 
from Tables 9; 10, only the median values of the three noise 
parameters were used. Therefore the error bars in noise Monte 
Carlo represent only the statistical variance in 101 simulated 
noise map realizations. If the uncertainty of the estimation of 
the three noise parameters was propagated to the noise MC, the 
error bars would be much larger.) 



Table 12. Summary of systematic effects uncertainties on maps" 
in ^Kcmb- 



30 GHz 



44 GHz 



70 GHz 



p-p rms p-p rms p-p rms 

Bias fluctuations 0.08 0.01 0.10 0.02 0.23 0.06 

Thermal fluctuations .. . 0.61 0.11 0.40 0.08 1.17 0.20 

1-Hz spikes 0.87 0.17 0.14 0.03 0.60 0.12 

Sidelobes pickup 18.95 4.53 1.92 0.57 6.39 1.91 

ADC non linearity 3.87 1.01 0.89 0.19 0.92 0.19 

Gain residuals 4.33 1.16 4.74 0.97 6.51 1.10 

Total 21.02 4.83 5.61 1.13 7.87 2.00 

il Calculated on a pixel size approximately equal to the average beam 
FWHM. 

12.5. Updated systematic effects assessment 

Analysis of know systematic effects in Planck LFI is reported in 
detail in Planck Collaboration III (2013). Here we include only 
the summary Table 12 where we list the r.m.s. and the difference 
between the 99% and the 1% quantiles in the pixel value distri- 
butions. For simplicity we referto this difference as peak-to-peak 
(p-p) difference, even though it neglects outliers. It nevertheless 
effectively approximates the peak-to-peak variation of the effect 
on the map. 

Our analysis (Planck Collaboration III 2013) has shown that 
systematic effect uncertainties, are at least two orders of mag- 
nitudes below the CMB temperature anisotropy power spectrum 
and are dominated by straylight pick-up from far sidelobes and 
imperfect photometric calibration. 



the results with the spectrum obtained from all 12 radiometers 
shown above. In Fig. 17 we show the difference between the 
horn-pair and the combined power spectra. Again, the error b ars 
have been estimated from the "FFP6" simulated dataset. A \ 2 
analysis of the residual shows that they are compatible with the 
null hypothesis, confirming the strong consistency of the esti- 
mates. 



13. Infrastructure Overview 

The computer cluster used for the maps production consists of 
two types of nodes: ten 64 bit nodes with two single-core CPUs 
and 16 GB of RAM, and ten 64 bit nodes containing two moth- 
erboards, each equipped with two six-cores CPUs and 72 GB 
of RAM. The total RAM available exceeds 1.5 TB, a sufficient 
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Fig. 15. Consistency between estimates of the angular power spectra at frequencies. From Left to Right: 70 vs 44 GHz; 70 vs 30 GHz; 
44 vs 30 GHz. Solid red lines are the best fit of the linear regressions, whose angular coefficients a are consistent with unity within 
the errors. 




Fig. 16. From Left to Right, the empirical distribution (estimated via "FFP6") of the s\, S2, S3 statistics (see text). The vertical line 
represents 70 GHz data. 



quantity to allow for the creation of all the maps until the end of 
the mission. The dual-motherboard nodes are connected through 
an infiniband network interface (40 Gbit), while an 1 Gbit inter- 
face is provided for the other connections. 

The hardware infrastructure includes a front-end (two quad- 
core CPUs, 8 GB of RAM), which is the access point for the 
users and hosts the pbs server, and a control machine running the 
LDAP authentication server and the DNS and DHCP services. 

The software used for the system management and synchro- 
nization includes tools as kickstart and puppet, while the 
parallelization of the computation is guaranteed by the torque 
resource manager and the maui scheduler. 

Data products are stored and organized into three different 
servers that host the Level 1, Level 2 and test databases (Fig. 1. 
For each of these databases there is an associated RAID 6 storage 
with up to 40 TB formatted with the JFS filesystem. 



14. Discussion and Conclusions 

We have described the pipeline used to process the Planck/LFl 
data, starting from Level 1 continuing the production of the tem- 
perature frequency maps based on the first 15.5 months of ob- 
servations. Furthermore, we have described the strategy for the 
verification of the quality of the products, which is largely based 
on null tests. Due to the complexity of the analysis process, three 
companion papers provide detail on specific critical aspects of 
the data analysis and products delivered: Planck Collaboration 
III (2013) reports the analysis on systematic effects and assess- 
ment of their impact; Planck Collaboration V (2013) describes 
in detail the photometric calibration (approaches and testing); 
Planck Collaboration IV (2013) outlines the determination of the 
LFI beam patterns and window functions. In addition, the Planck 



explanatory document Planck Collaboration ES (2013) provides 
a detailed description of all the products delivered in this release. 

The entire Level 1 pipeline was unchanged since the start of 
the mission (about four years ago) and has been running flaw- 
lessly and continuously, demonstrating the robustness of the de- 
sign and development approach. IN contrast, the Level 2 pipeline 
was largely restructured; see Zacchei et al. (201 1) for a descrip- 
tion of the initial pipeline. The major improvements involved 
new procedures for pointing reconstruction, detailed estimation 
of systematic effects and photometric calibration. These im- 
provements allowed us to obtain, as reported in Table 1, a final 
calibration uncertainty of the order of 0.6%, and also to propa- 
gate, using simulations, known systematic effects into the final 
product maps. The impact of the combination of all systematic 
effects has been evaluated to be at least two orders of magnitudes 
below the CMB temperature anisotropy power spectrum. 

Particular emphasis was given to null tests, which were rou- 
tinely applied to various subsets of the data, in order to quan- 
titatively assess the scientific quality of the LFI products. The 
null test procedure, described in Sec. 12, allowed us to detect 
and solve a number of problems in the implementation of our 
pipeline, which emerged during the processing period. In fact, 
the pipeline is still being optimized and more improvements 
are planned for the next data release. Future improvements of 
Level 2 will be aimed at obtaining high-quality polarization re- 
sults, which require control of spurious effects at sub-juK level. 
We are also concentrating our effort on a better beam charac- 
terization that takes in account second order effects - such as 
the bandpass response of each diode, or Galactic straylight (i.e., 
leakage through the far sidelobes) - during the calibration pro- 
cess. These and other refinements are being included in the LFI 
pipeline to meet the level of accuracy needed for a robust analy- 
sis of polarization data. 
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Fig. 17. Residuals between the auto power spectra of the horn 
pair maps and the power spectrum of the full 70 GHz frequency 
map. Error bars are derived from "FFP6" simulations. 
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